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1 Introduction 

The internal Diffusion Limited Aggregation (iDLA) model was first intro- 
duced by Diaconis and Fulton in |10| and gives a protocol for building a 
random set recursively. At each step, the first vertex visited outside the 
cluster by a random walk started at the origin is added to the cluster. 

The question of the limiting shape of this model is a well-studied one. 
Lawler, Bramson and Griffeath |21j were the first to identify the limiting 
shape of the model, in the case of symmetric random walks, as a Euclidean 
ball. Their result was later sharpened by Lawler [20] , and recently became 
of renewed interest with the simultaneous works of Asselah and Gaudillere 
(0, 0) and Jerison, Levine and Sheffield (P], (TTJ, [S] ). 

Still using simple random walks, Peres and Levine introduced in [23] 
a method for identifying limiting shapes other than a Euclidean ball, in 
particular the limiting shape of multiple source models. 

All these results are proved in the framework of simple random walks, 
and other random walks did not appear in the literature until Blachere's 
article [5j. In this paper, the iDLA model is studied for centered random 
walks, and convergence towards the ball of a specified norm is proved under 
moment conditions on the random walks. At the end of the paper, the case 
of drifted random walks is studied, and a limiting shape is found in the one- 
dimensional case. The author then conjectures the existence of a limiting 
shape for all dimensions. The initial idea that the limiting shape should 
be a level line of the Green's function, as it is the case in several types of 
groups (see [7], [6j) is disproved. 

In this paper, we present a limiting shape result for a simple class 
of drifted random walks. The limiting shape of the normalised cluster 
is characterised as a true heat ball because it gives rise to a mean-value 
property for caloric functions. The existence of such a bounded shape is an 
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open problem in PDE theory (see [H]), for which our convergence provides 
a proof inherited from the field of random walks. 
The main result in this paper is the following: 



Theorem 1.1 Let A n be the normalised drifted internal diffusion limited 
aggregation cluster. Then there exists a set D C ]R d_1 x M. + with the fol- 
lowing properties: 

1. Almost surely, A n converges towards D with respect to the Hausdorff 
distance. 

2. Let <p be a C°° function of time and space such that 

1 — p 



Then the following mean value property holds: 



(z,t)d(z,t) = \D\<f ) (0)- 



D 



3. The set D is bounded in time and space. 

The first point of the theorem states our convergence result for the 
normalised aggregate. The model admits a deterministic limiting shape 
whose properties are detailed in points 2 and 3. The second point states a 
mean value property for caloric functions, using the shape D, which justifies 
our use of the term true heat ball. The last point provides an answer to the 
as yet open problem of finding a such a mean value property on a bounded 
set. In a nutshell, the normalised iDLA cluster converges almost surely to 
a bounded true heat ball. Figure [T] is a simulation of this shape, obtained 
with 500 000 particles. 

Our proof follows the general idea of Levine and Peres, whose method 
for finding limiting shapes can be translated in our context. First, we 
introduce an equivalent of the divisible sandpile model (for a definition 
and convergence result of the original model, see [22J), which we call the 
unfair divisible sandpile, because one direction is privileged throughout the 
construction of the cluster. We study this model in details, and provide a 
limiting shape result. We prove convergence towards an abstract shape D, 
which for now is defined through a parabolic obstacle problem. 

Parabolic obstacle problems are relatively frequent in the literature, and 
occur in a number of varied fields. First and foremost, they are studied 
in the context of heat diffusion, and in particular in the case of two-phase 
transition equations, like the Stefan problem (see |23] and [9j for a dis- 
cussion of general parabolic free boundary problems). They also appear in 
finance, namely in problems related to the pricing of American put options: 
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Figure 1: Drifted iDLA aggregate with 500 000 particles and p = 0.2 

see for example [19] , [25], pj)]. In our case, we will use very strong results 
in parabolic potential theory (see for example [8], pt]) to ensure that the 
limiting shape D is smooth enough, and to characterise it as a true heat 
ball. 

The next step of our proof is to relate the drifted iDLA and unfair 
divisible sandpile, and prove that they almost surely share the same limiting 
shape D. To conclude, we use probabilistic arguments to give bounds on 
the drifted iDLA cluster, which are in turn used to bound D, thus proving 
the convergence of our normalised cluster towards a bounded true heat ball. 

2 Heuristics 

In this section, we will motivate the introduction of the divisible sandpile 
model for the study of the iDLA model. We consider the odometer function 
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introduced in |23j . that measures the total number of particles emitted 
from point x G Z, d throughout the construction of the cluster, counted with 
repetitions. 

A given point will then start with a certain number of particles (n if 
it is the origin, otherwise, in our situation), receive new particles during 
the construction of the cluster, and end with a new number of particles (1 
if it eventually lies inside the cluster, otherwise). 

Of all the particles that passed through a given point x, it seems natural 
that a proportion travelled to each of the reachable neighbors of x, and that 
proportion should in some sense be close to the transition probability from 
x to that particular neighbor. 

Assuming that the same holds for each of the neighbors of x, we get a 
tentative local equation: 

^2p(y, x)u(y) - u{x) = v{x) - a(x), 

y^x 

where cr(x) and v{x) are the initial and final amounts of mass at point x, 
respectively. 

While this reasoning is flawed because it assumes independence of cor- 
related quantities (among other problems), it motivates the introduction 
of a new model in which such a local relation holds. The unfair divisible 
sandpile model will play this role in our case. 



3 Unfair divisible sandpile 
3.1 Definitions and notations 

We introduce the unfair divisible sandpile model to be the drifted counter- 
part of the divisible sandpile model defined in [22J. Consider a continuous 
distribution of mass on Z d , with finite total mass and bounded support. A 
lattice site is full if it has mass at least 1. Any full site can topple by keep- 
ing mass 1 for itself, and distributing the excess mass among its neighbors. 
While in the classical divisible sandpile model, the mass is split equally 
among neighbors, in our model, it will be distributed proportionally to the 
step distribution of the drifted random walk S. 

At each time step, a full site is toppled. When we let the time go to 
infinity, if every full site is toppled infinitely often, the mass converges to 
a limiting distribution in which each site has mass less than 1. This is the 



object of lemma 3.1 



Just like the divisible sandpile model, our unfair version is abelian, in 
the following sense: while individual toppling do not commute, the limiting 
distribution of mass does not depend on the order of the topplings, provided 
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that they are done in an appropriate fashion, which is to say that each site 
that becomes full in the course of the topplings is then toppled infinitely 
often. 



We prove this abelian property in lemma |3.2| 

A crucial tool for the study of the divisible sandpile model, be it classical 
or unfair, is the odometer function. We define it as the mass emitted from 
x throughout the construction of the cluster: 

u(x) = total mass emitted from x. 

It is important to note that this quantity does not depend on the se- 



quence of topplings either, which will derive from our proof of lemma 3.2 



Since the mass emitted from a point is always distributed in the same fash- 
ion, we can remark that the mass received by x from his neighbors is the 
following: 



mass received in x = pu(x — e d ) H — — — ^— u(y) 

2 ( d 1 ) 

y~x,y-x±e d 



2(d 



Since we have defined u as the emitted mass, the difference between the 
received mass and u at point x will be equal to the difference between the 
initial and final amounts of mass at point x. Namely, 

u(x) - a(x) = pu(x - e d ) + * _ P ^ u {y)~ u { x )i 

y~x,y-x±e d 

where a and v are the initial and final amounts of mass. We split the 
quantity u so as to make sense of this equation: 

u(x) - a(x) = pu(x - e d ) + — - _ P u(y) - u(x) 

y^x,y—xXed 

= p (u(x - e d ) - u(x)) + * _ P ( u (y)- u ( x )) 

= —p(u(x) — u(x — e d )) + (1 —p)Au(x), (3.1) 

where A is the discrete Laplacian taken over the first d — 1 coordinates of 

Z d , which is to say that: Au(x) = J2 y ~ x , y - X ±e d ( u (v) ~ u ( x ))- 

We now define a discrete operator K, that sums up this operation, and 
will play an important role in our proofs: 

Kf{x) = (1 - p)Af(x) - p (/(*) - f(x - ei )) . 



We can now restate equation (3.1) in the following way: 

v(x) — a(x) = K,u(x) 
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Given the nature of /C, and since we will be inclined to consider its 
continuous counterpart, we will transform notations when we pass to the 
continuous: the first d — 1 will correspond to the space coordinates, while 
the last coordinate will correspond to the time coordinate t. The continuous 
counterpart of K, will be the following operator: 

1 — v df 
&f(x,t) = 2 (d- l) AfM ~ V Tt 

which is well known as the heat operator. Hence, we will henceforth call 
/C the discrete heat operator, and define a discrete caloric (respectively 
supercaloric) function as a function / : Z d — > R such that G r L d , Kf(z) = 
(respectively JCf(z) < 0). 

The sense in which ^ will be the continuous counterpart of JC will be 
a non-trivial point in our proof, and will give rise to the non-standard 
normalisation we introduce in section 13.3.11 

3.2 Convergence and abelian property 

We will now state the convergence result and abelian property for the 
unfair divisible sandpile model. Let us first define a toppling scheme. A 
toppling scheme T for a configuration Vq is an infinite sequence of indexes 
in Z d in which each site that is initially full or becomes full through the 
realisation of previous topplings appears infinitely often in the remainder 
of the sequence. We will denote by v£ the distribution of mass after the 
toppling of the first k sites listed in T (in the specified order), and u^(x) 
the mass emitted from x up to the k-th toppling. We omit superscripts 
when only one toppling scheme is involved. 

In addition, a toppling scheme will be called legal if it only tries to 
topple full sites. 

Lemma 3.1 Consider an initial configuration u with a finite total amount 
of mass M and bounded support S, and a legal toppling scheme T for this 
distribution. Then as k tends to infinity, and Vf. tend to limits u and 
v. The limiting configuration v is such that\/x e Z d , v(x) < 1. Moreover, 
these limits satisfy the following mass relation: 

v = u + JC(u). 

Proof: Let B be a ball centered at the origin big enough to contain all 
points within graph distance M of S. We remark that any point that if 
Vk{x) > at one point in the construction of the cluster, it means either 
that x was initially in the support of v$ or that it has received mass from 
a neighbor with mass greater that one. Since the same reasoning can 
be applied to this particular neighbor, it means that all points that have 
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positive mass must be within graph distance M of S. Hence, v k is supported 
in B. 

We introduce the weight function: 

W k = ^2 (^1 H h X d-1 + Xd ) ■ 

xez d 

Note that W can be negative, which does not present any problem as the 
key property is only that it increases through topplings. 

When toppling site x at time k, the mass at point x has been modified 
by an amount a k (x) = — v k (x) which has been transferred to the 

neighbors of x according to the toppling rule. Hence, 

w^-w^x) = a k{x) ^~v yi-xl + --- + y 2 d-i-x 2 d -i 

y~x,y-x±e d 

+pa k (x) (x d + l- X d ) 
= a k (x) (3.2) 

Since u k is the sum up to k of all the relevant eti(x), we get the following 
relation on weights: 

W k = W + 

x& d 

For every x, u k (x) is a bounded increasing sequence, so it converges to a 
value u(x). 

At all finite times k, the relation 

u k (x) = u (x) + JCu k (x) 

holds, so that the convergence of v k is a consequence of that of u k . More- 
over, its limit v is such that v = is + /Cm. 

Now a point x is either never toppled, in which case we have v k {x) < 1 
for all values of k, or it is toppled infinitely often, but then each time a 
toppling occurs at point x, u(x) < 1. In both cases, the limit v{x) has to 
verify u(x) < 1, which concludes the proof. 

We next prove the abelian property of the unfair divisible sandpile 
model. 

Lemma 3.2 Consider an initial configuration i/ , and two legal toppling 
schemes S and T. Then the corresponding final configurations v s and v T 
are equal, as are the final odometer functions u s and u T . 

Proof: Let x k be the point at which the fc-th toppling occurs in the scheme 
S. We will prove by induction on k that u T (x k ) > u k (x k )- 

Suppose that the property holds for all i < k. Then for x ^ x k , either 
x is not toppled before time k in scheme S, and u k (x) = 0, or it is, in 
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which case we consider the last index i for which toppling occurs at point 
x. Then u T (x) > uf(x) = u k (x). In both cases, u T (x) > u k (x). 
Since S is legal, and T is a toppling scheme, 

V T (x k ) < 1 < I/fc(Xfc). 

Hence, equations (3.1) and ( 3.2[ ) ensure that: 

K,u T (xk) < K.u (xk). 

We write the mass relations for point x k . 

u T {x k ) - fCu T (x k ) = u (x k ) - v T {x k ) 
u s k {x k ) - K,u s k {x h ) = v Q {x k ) - vl(x k ). 

Taking the difference yields: 

u T (x k )-u s k (x k ) > 2 ^~ p ( uT (y)- u i(y)) 

y~x k ,y-x k ±e d 

+p (u T (x k - e d ) - uf(x k - e d )) 

> 0. 

The right hand side is indeed positive since it only involves differences for 
points different from x k . 

Hence, we have proved by induction that for all k, u T (x k ) > u k (x k ). It 
follows that U2 > U\. A symmetric argument shows that U\ > u 2 , which 
concludes the proof. 



3.3 Limiting shape 
3.3.1 Normalisation 

Since the normalisation we will use is non-standard, we introduce a specific 
normalisation function <f> n . The fact that this normalisation is not the same 
in the direction of the drift as in the other directions plays a very important 
part throughout the paper. It is defined as follows, for functions defined 
on Z d with values in ffi. : 

V/ : Z d -»■ R, (j) n {f) : R^ 1 xl^R 

(f> n {f)(x,t) = f ((Ln^(s)J) i6{w _i } , [n^t\J . 

With a slight abuse of notation, we extend the definition of n to subsets 
A of Z d : 

<Pn{A) = { (x, t) G R^ 1 x R, (( [nsri ( X )A ) ie{1 ,..., d _ 1} , [n^t\ ) G a} . 
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3.3.2 The heat equation 



In the following sections, we will be required to handle the heat equation, 
and in particular to define what a supercaloric function is in the continuous 
setting. We will follow exactly the definition of Evans, [IT]. First, we define 
the fundamental solution <&(x,t) of the heat equation R(f) = as follows: 



$0,t) 

where we define 




t > 0, 
t < 0, 



2(1 -p)' 

Note that $ is singular at the origin, and solves = away from the 

origin. 

For fixed x G E d_1 , t G E, r > 0, we define: 

ul i - - ■ ■■■ ' 



E {x ,t-r) = { (y, s) G M x M, s < t, $(s — y,t — s) > 



d-l 



Note that the boundary of is a level set of $(x — y, t — s). Point (x, t) 
is at the top, center end of E. 

Let / be a measurable function; we define the operator 2l r as follows: 



%r(f)(x } t) = I f(y,s) j V V s - dyds. 



- l (x,t;r) 

This operator is a sort of mean value operator; indeed, if / is a smooth 
solution of the heat equation, 2t r can be used to compute the value at a 
given point: 

u(x, t) = %l r (u)(x, t). 

Note that the right hand side involves only u(y, s) for times s <t, which 
is reasonable, because the value at a given time should not depend upon 
future values of u. 

We now use this operator to define supercaloric functions, as follows. 
Let / be a measurable function. We say that / is supercaloric if it verifies 
the following property: 

V(x,t) G M^ 1 xl,Vr> 0, f(x,t) > %.(f)(x,t). 

For a function that has C 2,1 regularity, this is equivalent to the intu- 
itive result: 

< 0. 



Note that equation (3.3.2) is only a weak substitute for the mean- value 
property of harmonic functions; instead of simply integrating the function 
on a given shape, a kernel is used. The problem of the existence of a gen- 
eralized mean value property for caloric functions, which is a consequence 
of Theorem is discussed in 
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Figure 2: The obstacle function j(x,t). 
3.3.3 Main result 

The following theorem provides a result for the limiting shape of the unfair 
divisible sandpile aggregate under the proper normalisation. 

We run the unfair divisible sandpile on 7L d with the following initial 
mass configuration: 

o n {x) = nS {x), 

and we intuitively define the unfair divisible sandpile aggregate as the 
set of points in Z d that have positive final mass: 

D n = {xeZ d ,u n (x) >0}, 

where v n is the final mass configuration corresponding to the inital mass 
distribution /i n . 

Recall that (3 = 2a-p) > an< ^ define the following quantities: 

7 0M) = t -\\x\\*-±(jLy exp(-/3&), 

s(x,t) = mf{h(x,t)\h(x,t) is supercaloric and h > 7}, 
u(x,t) = s(x, t) — j(x, t) 

D = {(i 1 t)6i"xR;,!i(i 1 i)>0} = {s-7>0}. 

The function 7 is plotted in Figure [2] when d = 2. 

Then, the following theorem states the convergence of the normalised 
unfair divisible sandpile aggregate on compacts. 

Theorem 3.3 For every compact K ofM. x ]R^_ ; the restriction to K of 
the normalised unfair divisible sandpile aggregate <p n (D„) DK converges to 
D n K with respect to the Hausdorff distance. 
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We start by proving that the odometer function for the unfair divisible 
sandpile with starting configuration k at the origin is given by the differ- 
ence between an obstacle function 7 n (for which we fix JCj n ) and its least 
super caloric majorant. 

Lemma 3.4 Start with mass n at the origin and choose 7 n such that 
/C7 n (x,t) = 1 if (x,t) 7^ (0,0), and /C7 n (0,0) — 1 — n. Then the odometer 
function u n is given by u n = s n — j n , where s n is the least supercaloric 
majorant of 

Proof: We want to prove the following: 

u n + 7n = inf{/(x,t)|/C/<0and/> 7n } (3.3) 

Let us first remark that u n + 7 n is indeed a supercaloric function such 
that u n + 7„ > 7„. 

We recall the mass equation (3.1 ) as stated in terms of the discrete heat 
operator: 

v n {x) - a n (x) = fCu n (x). 

Recall that the initial distribution of mass is n at the origin, while 
the final distribution has values that are everywhere less than 1. Hence, 
u n + In is supercaloric everywhere, and since u n > by definition, it is also 
a majorant of j n . 

Let us now prove that is it the smallest of these functions. Let / be a 
supercaloric majorant of j n . Then: 

£ (/ - 7n - «n) = JCf -K, (u n + 7) 

On D n , K (u n + 7„) = 0, so that /C (/ — 7„ — u n ) < 0. Outside D n , 
u n = 0, so that we have / — 7 n — u n = f — 7„ > 0. 

Hence, / — 7„ — u n is positive outside D n and supercaloric inside D n . 
The minimum principle yields that / — 7 n — u n > everywhere, and / is 
indeed always greater than u n + j n , which concludes the proof. 

Remark: Lemma |3.4 gives us a way to find the odometer function for 
the unfair divisible sandpile model, provided that we have a function that 
satisfies the following conditions: 

/C7„(x,t) = 1 if (x,t) 7^ (0,0), and 

/C 7n (0,0) = l-Ti. (3.4) 

We introduce the random walk S on Z d , with the following law: 

W(S(t+l)-S(t) = ±ei) = 2{d P \) fori = l..-d-l, and 
¥{S(t+l)-S(t) = e d ) = p. 

(3.5) 
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Our function 7„ is easily written in terms of the random walk S. Define 
g(x, y) as the expected number of times that S, started at x, visits y. Then 
consider the function: 



Jn(z) = Zd-^zl - ng(0, z). 



i=i 



(3.6) 



It is easy to check that this function verifies the conditions (3.4). 

In order to be able to use this function, we will need to evaluate the nor- 
malisation of Green's function, which is the object of the following lemma: 

Lemma 3.5 For z G 7j d , Define g n (0,z) = ng(0,z). Then the normalised 
version of g n , namely n d + 1 (fi n (g n ), converges uniformly on compacts of 



l d 1 x ]R*_ towards the function 



G(x,t) = - 1 - ( P-) 2 expf-/3 lklh 



p \nt J \ t 



To prove this lemma, our first step is to couple our random walk S with 
the random walk S, whose steps are the sum of those of S between two 
jumps in the direction of the drift. In other terms, the law of the increments 
of 5* is that of the last position of simple random walk killed at a geometric 
time of parameter p. 

We calculate the covariance matrix of the walk S, so that we can esti- 
mate its position at a given time using the local central limit theorem (see 
[26], P9). 

Let us denote Pfc(0, z) — P = z^j. Then applying the local central 
limit theorem, we get that for all (x,t) G ]R d_1 x W + , 



d-1 



P\\x\ 
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na+r (27rt) ^ P [n ^ ti (0, [n^x\ J (2/3)^ exp f 

and the convergence is uniform on compacts of x M.* + . Moreover, the 
difference between these quantities is of order ^~ 3 + T ^ ) see f° r instance 
1131. 



The probability we estimated here is that of the event that the first 

2 1 

site visited on the L^-^+^J-th layer in the drift direction by S is L^ 3 ^^- 
Hence, 



g(0,(n^x\,[n^t\^ = £ P ^ (o, [n^x\ + g{(z, 0), 0) 

(3.7) 
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Equation (3.7) states that in order to get the actual number of visits 
to |_rz rf+i , we have to consider the sum of the hitting probabilities of 
[n^+ixj + z, where z G Z d_1 , multiplied the number of visits to of the 
drifted random walk started at point (-2,0). The hitting probabilities of 
[n^TTa;] + z with z small enough to be in range of are asymptotically 
close to that of [ri^+^xj , so that we only have to evaluate the following 
sum: 

s(M),o) = 

using the fact that the origin and (z, 0) both have coordinate in the drift 
direction. The summation yields exactly the expected time spent by the 
random walk S on one given layer. 



3.3.4 Convergence of the obstacle function 

2 

Lemma 3.6 The normalised obstacle n~^+ x (j) n { n f r} ) converges uniformly on 
compacts K ofM d ~ x x R towards 7. 

The proof of this lemma is straightforward since the first two terms in 7 n 
converge uniformly to those of 7 once normalised. The last term is the 
normalisation of Green's function, the convergence of which was proved in 
lemma 13.51 

2 

Hence, our normalised obstacle function n~ 3+ 1 0„ (7„) converges uni- 
formly on compacts toward the function: 

7(*,t)=t-|M| 2 -^) 2 exp(-/*ML). 

3.3.5 Parabolic Potential Theory 

Before we go on, we will need a few lemmas of parabolic potential theory. 
Since n~d+i(j) n (g n ) converges uniformly on compacts of ]R d_1 x M.* + towards 
the function 

G <°^ = -K£>" exp (-^r)' 

we want to use this convergence to find a candidate to be the inverse of the 
heat operator, which we define in both the continuous and discrete setting. 
Let us define, for zGZ^ and / a measurable function defined on M?- 1 x R* 
with compact support, 

G n {f){z) = n 9( z ,v)f (n~^y,n-£ir) , 

(y,r)eZ d ,r>0 
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and for x G R^ 1 x R* + , 

G(f)(x) = [ G(x,(y,r))f(y,r)dydr. 

J (y,r)£R d - 1 xR* + 

The uniform convergence of the normalisation of g n towards g, together 
with our estimate, will ensure the convergence of the normalisation of G n (f) 
toward G(f). This is the object of the following lemma: 

Lemma 3.7 Let f be a bounded function defined on W* -1 x R*^ with com- 
pact support, then 

|n-^ n (G n (/))-G(/)|->0, 

uniformly on compacts ofM. d ~ 1 x R+. 

Suppose that / is bounded by M and supported on the compact K. 
We define, for (y,r) G K such that n (y,r) G 7L d , the following error term: 



<(y,r) = 



_ 2 

n d +! 



<f> n {g) (x, (y,r)) - / G(x,z)dz 

J z 



where the integral is taken over points z G (y, r) + [0, n d +! ) d 1 x [0, n d + 1 ). 
Then we have: 

n-^UGn(f))(x)-G(f)(x) < J2 Ma n(v> r ) ( 3 - 8 ) 



where the sum is taken over points (y,r) G K such that <j) n (y,r) G r L d . 

__ 2 

Since we have a uniform bound on the difference n d +! (f> n (g) (x, (y, r)) — 
G(x, (y,r)) that has finite integral on K, and using our knowledge of the 
derivatives of G to bound G (x, (y, r)) — G (x, z) when z is close to (y, r), 
we can bound this sum uniformly on K, which concludes the proof. 

The idea behind the introduction of G(f) is that G is in some sense the 
inverse of the heat operator. In fact, we do not need such a strong result, 
and we are content with the following property: 

Lemma 3.8 Let f be a positive measurable function defined on R 4 - 1 x R* 
with compact support, then G(f) is supercaloric. 

We omit the proof as it is straightforward and only relies on commuting 
integral signs and using the supercaloric property of G(x,y) as a function 
of x. 

We are now ready to prove the convergence of the odometer functions 
in our model. 
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3.3.6 Convergence of the odometer function 

Let us consider s n , the least supercaloric majorant of j n . The following 
lemma states a convergence result of the normalised version of s n . 

Lemma 3.9 The normalised version of the least supercaloric majorant of 

2 

7 n; namely n~ 3 + T ri (s n ) ; converges uniformly to s, the least supercaloric 
majorant of~f, on compacts ofR d ~ l x R* + . 

We begin our proof by pointing out that our definition of the least 
supercaloric majorant depends on the context: while s n is defined in the 
discrete setting, s is defined in the continuous. 

Let K be a compact of IR d_1 x R* + , and Kj an increasing sequence of 
compacts of R^ 1 x R* such that UjKj = R^ 1 x R* + , and Kj C K° j+1 . We 
also define their discrete counterparts K™ as the set of points (x, t) such 

that (n~^x, n~^tj G Kj. 

We define the following quantities: 

s£ = inf {f{x, t) \Kf <0 globally, and/ > 7n on K?}, 

s Kj = mi{h(x, t)\h(x, t) is globally supercaloric and h > 7 on Kj}. 

Remark that both functions are increasing in j, so it is a consequence 
of Dini's theorem that s Kj and 0„(sn J ) converge uniformly on K towards 
s and (f> n (s n ), respectively, as j tends to infinity. 

We now write that: 

\MSn) -S\< \4> n {s n ) - <t> n {s%)\ + \<t>n{s%) ~ | + - s\ 

Set j big enough so that K C -Kj, and \(f) n (s n ) — n (sn" J )| + — s\ < e 
on i^. It now remains to show that ifinisn 3 ) converges uniformly on K 
towards s Kj . 

We will now proceed to smooth the function s Kj in order to show that 
it is close to a function s K: > such that JCs Kj is small. 

We first remark that s Kj is continuous on K. Indeed, as it is an infimum 
of continuous functions, it is also upper semi-continuous. Moreover, if we 
define a; (7, r) as the continuity modulus of 7 on Kj +i , we have, on Kj, 

% r (s Kj ) > 2t r (7) > 7-w(7,r), (3.9) 

so that 2t r (s Kj ) + a; (7, r) is continuous, supercaloric, and lies above 7 on 
Kj. Hence, 

2t r (s^) < s K > < %. (s K i) + w( 7 , r). 

Since 7 is uniformly continuous on Xj+i, 2t r (s^ J ) — > s Xj as r — )■ 0. 
Since is supercaloric, this is an increasing limit. As an increasing limit 
of continuous functions, s Kj is also lower semi-continuous. 
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We now define, like in [TT], Appendix C, for x G K, A > 0, 



where rj is defined as follows: 

{x) = J Cexp , N<1 

I 0, |ar| > 1, 

with C such that J Rd r/=l and A is taken sufficiently small to ensure the 
definition of the values of s K * involved. Then s K] is still supercaloric, and 
if A is small enough, it verifies \s Kj — s K: >\ < e on K. 
Let us define the following discrete function: 

/ 1 2 \ 5 

(yk,l) E Z^ 1 x Z, q n {x,t) = s K] (n _ 3+Tx,n"3+Ttj - | x |2_ 

Since s Kj is smooth and supercaloric, it verifies the partial differential in- 
equation: 

< 0. 

Moreover, since s Kj is smooth, Taylor's formula yields that, on K, 



. _3_ 

< An d + 1 . 



where A is chosen to be the maximum of the norms of third derivatives in 
space dimensions, and second derivative in time of s j . 

Hence, the smoothness and supercaloric property of s Kj ensure that, 
for n large enough: 

JCq n (x,t) < o. 

When n is large enough, q n is close to s j , which is in turn close to s Kj , 

2 

which is greater than 7, itself close to ri _2 + T n (7 n ) on Kj. To sum up, the 
following inequalities hold on Kf 



<Pn(q n ) > s Kj - e > s K > - 2e > 7 - 2e > n <w 4> n (ln) ~ 3e. 

2 

Hence q n + 3e is a supercaloric majorant of n~"^+ i <p n { n f n ) on Kj, so 

Ms*') < MQn) + 3e < ~s K > + 3e < s K i + 4e. 
We will now prove the converse inequality. 
Define the following functions: 

G(hj)(x) = / G(x,(y,r))hJ([nW : y\,[n'^r\)dydr. 
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We will start by proving that on K^ +1 , 



To see this, a discrete reasoning is necessary. Indeed, let us look at one 
particular point z, and suppose K,s n J ' (z) < — 1. Consider the function / 
that coincides with Sn 3 on every point but z, and adjust the value f(z) 
such that Kf{z) = — 1. Remark that this implies f(z) < Sn 3 (z). We will 
now prove that / is a supercaloric majorant of 7„. First, it is supercaloric 
at point z because JCf(z) = — 1 by definition, and on other points, since 
f(z) < Sn J (z) implies that for y ^ z, Kf(y) < K,Sn 3 {z). Note that, as an 
infimum of discrete supercaloric functions, Sn 3 (z) is itself supercaloric, so 
that for y z, JCf(y) < K,s n 3 {z) < 0. The function / is also a majorant of 
j n , which is true at points y ^ z by definition. At point z, we know that 
JCf = —1 = ICjn, so that the minimum principle applied locally between 
z and its neighbors to the function / — 7„ guaranties f(z) — J n (z) > 0. 
Hence, / is a supercaloric majorant of 7„ on K™ +1 , so that s n 3 < /, which 



is a contradiction. Hence, n K™, 



\Ks 



< 1. 



On the other hand, since — G n inverses /C exactly, n d + 1 s'n 3 
on all non-boundary points of K™ +1 . 

We want to argue that G(h^) (x) is a good approximation of n~ 3+r s 
We separate the error as follows, for x G Kj", 



G n ih n A 



where 



G{h n A{x) -rC^s%{x 



<A + B 



A < 



\ h 7(y^ r )\<(y^ 

2_ 

B < B n d +! 



(y,r),(z,t) 



Sn 3 ' {y, r)a£(z, t) - s% (z, t)a x n {y, r) 



where the second sum is taken over points (y, r) ~ (z, t) such that (y, t) G 
but (z,t) ^ -^j+i) an d B Q is a suitable constant depending only on 
the dimension. The term B estimates the part of the error that stems from 

in the definition of W. To bound A, we use 



using Sn'lic?, , rather than s n 



equation (3.3.6): 



A< <^ r ) 



We then find the same sum as in the proof of lemma 3.7 
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Bounding B is a little trickier. We proceed in the following way: 

K J (y,r)<(M) -s%(z,t)a*(y,r)\ < \s^(y,r)\ |<(z, t) - <(y, r)\ 

+ \<(y,r)\ \s^(y,r)-s^(z,t)\ 

We will rely on three arguments: first, remark that both (y, r) and (z, t) 
are safely away from x, because of our condition Kj C with (y, r) and 
neighbors, a short computation shows that: 



a%{z,t)-a£{y,r) = o (n 



Secondly, it follows from the definition of s n 3 and the convergence of 

2 

ra~^0 n (7„) that on K™, 



&{z)<Cn&, 

where C is a constant depending on Kj. 

Last, we use the fact that 7 is uniformly continuous on Kj to argue 

that the increments of 7 n between two neighboring points are at most of 

2 

order en d + 1 for n large enough, and this property is in turn transferred to 
s n 3 (since s n 3 {■ + e$) + e is a supercaloric majorant of 7 n ). 

Hence, our terms are bounded in a satisfactory way, so that A and B are 
asymptotically close. Hence, on Kj, G(h n ) is uniformly close to rT^Sn 3 , 

2_ 

which is bigger than n d+1 7 n , which is uniformly close to 7. To sum up, 
G(h n ) > n~J^is% J - e > n~tt"f n - e > 7 - 2e. 



on i^j, 



Moreover, lemma 3.8 states that G(h n ) is supercaloric everywhere, so 



that G(h n ) + 2e > s and we can conclude that Sn J + 3e > s Kj . 

The convergence of s n 3 towards s Kj is thus uniform on K. Recall that 
the convergences of s Kj towards s and s n 3 towards s n 3 are guaranteed to 
be uniform by Dini's theorem, which concludes the proof. 

2 

Lemma 3.10 The normalised odometer function n~ ^+1 (f) n (u n ) converges 

towards u = s — 7 uniformly on compacts o/R d_1 x R^_. 

2 

Moreover, the function u and the family n d+1 4>n(u n ) have bounded (re- 
spectively uniformly bounded) integrals over any compact of x M + . 



The first statement is a consequence of the convergence properties of s n 
and 7„, while the second relies on the fact that t — \ \x\\ 2 is a supercaloric 
majorant of 7 (resp. j n ). Therefore, for x G Z d , 

s n {x) - 7„(x) < ng(0,x), 
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and for (x,t) G R^ 1 x R* + , 

A quick computation finishes the proof. 
3.3.7 Convergence of Domains 

Lemma 3.11 Let K be a compact ofR^ 1 x M* + . Then, on K, the nor- 
malised unfair divisible sandpile cluster <p n (D n ) converges to D with respect 
to the Hausdorff distance. 

Set e > 0, and define D e as the inside e-neighbourhood of D, that is to 

say, 

D e = {xe D,B x , e C D}. 

Since D is the non-coincidence set for the caloric obstacle problem with 
obstacle 7, s— 7 is strictly positive on D e . Since s— 7 is uniformly continuous 
on K fl D e , let us define 

j3 = min_ (s(z) — 7(2)) > 0. 

Then, for n large enough, the uniform convergences the normalised 
obstacle and majorant guarantee that I7 — n _3 + T ri (7 n )| < |, and \s — 
n _3 + I n (s n )| < | on K fl l). This yields the following : 

n"^ T 0„(s„) - n~^T0„(7 n ) >s-7-/3>0on^nD e , 

so that K C\ D e <Z K C\ n (_D( n )) for all but finitely many n. 

Let D e be the outside e-neighbourhood of D. Let (x ra ,t n ) G Dr n \, and 

suppose that (n~d+T;r n , n*+it n ) converges to (xo,t ) G M.^ 1 x It is 

sufficient to prove that (rc ,to) G 

Define the discrete cylinder segment C = B (x n) n^ 1 1 j x [ [_t n — '^*+ T , \ t n \ ] ■ 

On C fl £>(„), define the following function : 

w n (x, t) = u n (x, t) - (\x - x n \ 2 - (t - t n )) 

Then JC(w n ) = Ku n — 1 > 0, so that — w n is discrete supercaloric. 

Hence, w n satisfies the parabolic maximum principle and realises its 
maximum on the parabolic boundary of C fl -D( n ). 

Since w n (x ,t ) = u n (xo,t ) > 0, and u n = on <9-D(„), the maximum 
cannot be taken on dD( n y 
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Hence w n takes its maximum on the parabolic boundary of C, d p C. This 
set can be decomposed as follows: 



d p C = Ai U A 2 , 



where 



2 

Ai = B (x n , 1) x { [tn - e —\ }, and 
A 2 = <9B(x n ,n^ 



_ | x 

2 ' 



Now, on A 1; we have: 



2 e 



w n (x,t) < u n (x,t) - n d +! — , 



and, similarly, on A 2 , 



2_ e 



w n (x,t) < u n (x,t) - n d +i 
Let (y n , r n ) be the point where w n takes its maximum. Then, 

< u n (y n ,r n ) 

This gives the following: 



2 £ 

«n(j/n, r n) > U n (x n , t n ) + — . 

Let us extract a subsequence of (y n ,r n ) that converges to a point (y, r). 
Since n'*+ T n (-u n ) converges to w uniformly on compacts, the limit u(y,r) 
is strictly positive, so that (y,r) e -D. 

However, point (y,r) is within distance e of (x ,to)) so we can conclude 
that (x ,t ) G D € , and D (n) C £> e . 



3.4 Regularity of D 

3.4.1 Differential Equation approach 

So far, we have considered u only as the limit of the discrete odometer, but 
it will be useful to see that it is also the solution of the continuous version 
of the equation that defines u n . This is the object of the following lemma: 
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Theorem 3.12 The normalised limit u of the odometers solves the follow- 
ing partial differential equation in the weak (distributional) sense: 



1 — p du 

W^T) P di = lu>0 

where 5q is the Dirac measure at the origin. 



-6n 



(3.10) 



Consider a discrete function rj that has compact support on Z d . Since 
u n solves the discrete equation K,u n = l Un >o — i^Sq, we have: 

Y,-niz)K.u n {z) = ^(^l^o- n77(0,0) (3.11) 

z&L d z£Z d 

When changing the variables so as to report the operation on 77, only the 
sign of the drift coordinate is modified, so that we define the discrete op- 
erator K, as the operator /C with reversed time: 



P 



2(d 



E 



(m-m)+p(m-f(z-e d )). 



\\y-z\\=i,(y-z)±e d 

Equation ( |3.11[ ) can now be written in the following format: 

J2un(z)iCri(z) = X>(z)l u „>o- ^7/(0,0). (3.12) 

z& d z& d 

Let us now consider a test function h : IR d_1 x R — > R that is C°° and 
has compact support. Define its discrete counterpart h n : Z d — > R as: 



h n (z 1 ,...,z d ) = hl(n d +iZi) ,n d +^z d ). 

\\ / ie{l,...,d-l} J 



Then applying equation ( |3.12 ) yields: 

1 x ~ 1 . 

- u n (z)K,h n (z) = - Y] h n (z)l Un>0 - h(0, 0). 

zez d z& d 

- (rT&Unizj) (n&Khnizj) = - J2 h n( z ) - M0,0). ( 3 - 13 ) 



zeD n 



Theorem 3.3 proves that <p n (D n ) converges to D on compacts K C 



and lemma 



3.10 



proves that n d+1 (J) n (u n ) converges to u uniformly on 



K. Hence, letting n go to infinity in (3.13), the regularity of h and its 



derivatives, together with the bounded integral property stated in lemma 



3.10 ensures that we can apply the dominated convergence theorem, which 
yields: 

(1 — p) A , dh\ 



u 



2(d-l) Ah + P d-t) 



h-h(0,0), 



D 
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and equation (3.10) holds in the distributional sense. 



This characterisation of u as a solution to a parabolic free boundary 
problem yields numerous analytical results on the function u and the set 
D = {u > 0}. The following proposition states the most important ones 
for our problem: 

Proposition 3.13 Let <fi be a C°° function of time and space such that 

1 — p A , dd> 

Then the following mean value property holds: 

<p(z,t)d(z,t) = \D\(p(0). 

D 

Moreover, the limit u of the odometer is continuous in both time and space 
coordinates, and the set dD has (d- dimensional) Lebesgue measure 0. 

The first part of the proposition is a direct consequence our characteri- 



sation of u as a solution to the equation (3.10) in the distributional sense. 
The second part of the property states regularity results on such a solu- 
tion. For the proof of this powerful result in parabolic potential theory, we 
refer the reader to j8j. For an in-depth study of parabolic free boundary 
problems and their regularity, see the detailed book by Friedman \V2\ . 

4 Drifted iDLA 

In this section, we prove the almost sure convergence of the normalised 
drifted iDLA cluster towards D, the limiting shape of the unfair divisible 
sandpile model. For reasons of simplicity and lisibility, we will start by 
using a random walk whose restriction to the direction of the drift is an 
increasing process. 

4.1 The model 

Let (S^'gN a sequence of independent random walks on Z d , with the fol- 
lowing law: 

F(S j (t + l)-S j (t) = ±ei) = }~ P , for z = l---d- 1, and 
v 2 (a — 1) 



F(S j (t + l)-S*{t) = e d ) = p. 



(4.1) 
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The iDLA cluster is built recursively using the random walk S j in the 
following fashion. We start with an empty set. At step j, suppose that we 
have built A(j — 1). We launch the random walk at the origin, which we 
let evolve until it exits A(j — 1). The first site visited outside the cluster 
is then added to the cluster. 

Formally, let us define the cluster A{n) and stopping times (o~k)km 
recursively in the following way: 

vi = 0, 

A(l) = {0} = {S 1 ^)}, and 

Vj > 1, uj = infO > : S j (t) £ A(j - 1)}, 

A(j) = A(j - 1) U {&(*;)}. 

A few simple facts can be stated about the cluster. It is a random 
increasing set that contains the origin. We have %A (j) = j, and given the 
form of the law of the random walks, the intersection of the cluster with 
jd-i x is empty. 

4.2 Limiting shape 

We first prove internal convergence of our cluster intersected with compacts 
of ]R d_1 x M* + . We then use this convergence, together with the fact that 
it yields control over a proportion of the particles arbitrarily close to 1, to 
prove external convergence as well. 

Theorem 4.1 (Internal convergence towards D, on compacts) Let A{n) 
be the drifted iDLA cluster, and K a compact subset ofW 1 ^ 1 x W + . 
Then, for every e > ; 

for all but finitely many n, almost surely. 

Recall that we build the cluster using the family of independant random 
walks (S l ) ie {i... n }. For z G Z d , we efine the following stopping times: 

tI = inf {t > 0, S*(t) = z} 

= infO>0,^(*),^ £> n } 
z/ = M{t > 0,5*^) i A(i- 1)}, 
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We introduce the following counting random variables : 

n 

£n(z) = ^2 1 u l <Tl<rl n 
i=l 
n 

M n {z) = Y.U<rh n 
1=1 
n 

i=i 

Since measures the number of particles that pass through z before 
either adding to the cluster or leaving D n , if J\f n (z) > 0, then point z is in 
the cluster A(n). 

Remark that M n {z) = M. n {z) — C n (z). 

The strategy of the proof will be to study these random variables. After 
showing that their expected value is different enough, and that they are 
close enough to their respective expected values, we will deduce that their 
difference is suitably big as n tends to oo. 

First, we use the Markov property of our random walk to argue that 
summing over all walks that add to the cluster before exiting D n can be 
re-written as a sum over points of D n . Let us construct the following family 
of random walks S y , for y G D n as follows. If there is an index % e {1 • • • n} 
such that y = S l {i>i), then we define S y (t) = S l (t + v^). Remark that the 
aggregation rule guaranties that there is at most one such index. If there 
is no such index, then we define S y as a random walk started at y, with 
the same increments as S* 1 for instance, and independent of all the other 
random walks. 

We define the corresponding stopping times: 



T y = M{t>0,S v (t) = z} 
rl n = in£{t>0,S»(t),tD n } 

(4.2) 



Then we have the following inequality: 

an 



i=l 

yeD n 



which holds because every term in the first sum can be found once in the 
second sum, with possible additional terms. Remark that C n (z) is now a 
sum of independent indicator random variables. 
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Let us compute the expectations of these variables: 
E(M n (z)) = n ^fA 

9n,D n {Z,z) 

where g n ,D n is the Green function of a random walk S stopped when it exits 
D n . 

Define T^n,D n the time at which S^, a random walk with opposite drift, 
starting at z, will exit D n . Then we have: 

c E{T ^\ < E (C n (z)) < C E{T ^\ 

9n,D n {Z,z) V J 9n,D n {z,Z) 

Define f n ,D n { z ) = 9n,D n (z,z)E (M n (z) - C n (z)^j . Then, 

fn,D n (z) = 9n^Z,z)(it ¥ ( Tl z< T k)-J2 F ( T z <T k) 

\i=l yeD n 

= g n , Dn (z,z) (*>(v)" - l)P(r» < r£J. 
yeD n 

= ^(^(y)n-i)gn,D n hji z ) 

yeD n 

Hence, f n ,D„ is a solution to the discrete parabolic PDE : 

fcfn,D n ( z ) = l-nSo(z) VzeD n , 

fn,D n (z) = *Z$D n . 

Recall that u n , the odometer function for the unfair divisible sandpile, 
satisfies the same equation, on the interior of D n . This means that f n ,D n — 
discrete caloric function on the interior of D n . Moreover, for z on 
the interior boundary of D n , since z has at least one neighbor who never 
toppled, it means that the total mass emitted from z towards this neighbor 
is less than 1, which means we have: 

u n {z) < 2d 

Hence, for all z G Z d , 

fn,D n {z) > Un(z) - 2d 

Let us define the following minimal value of u: 

8 = inf u(x). 

xeD t r\K 
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Since u is continuous and K is a compact, the infimum is actually a 
minimum and, we have (5 > 0. 

2 

Since n d + J 4> n (u n ) converges uniformly to u on K, we can choose n large 
enough so that for all points z G Z d such that: 



i 



ie{i-d-i} 



we have, for n large enough, 





u n (z) > n« 



This means the following relation on the expected values of C n (z) and 
M n (z): 

2 

/ \ 1 6n d + 1 

E (M n (z) - jC n (z)) > - P (4.3) 



Using the bounds in equation (4.2), we get: 



EiMJz)) > fl + ^ n ' +1 . ^)e(cJz)) (4.4) 

I 2C^ Ai (^,z)E(r- >D J / v nl V 



i 



Since n d+1 I , n d+1 ^ 1 G -D e D if, we have for a suitably 

\\ /ts{l-d-l} / 

large constant C", uniformly in 2, 



Hence, equation (4.4) can be written as 



E(M»(z)) > (1 + k)e(A»(2)), 

where k is a strictly positive constant. 

It follows that we can choose A > such that 

(1 + A)E (£ n (z)) < (1 - A)E (M n (z)) , 

and apply the following large deviation principle, which holds for sums of 
independent indicator variables : 

P (£ > (1 + A)E (£)) < 2e- c ^ 
P (M > (1 - A)E (X)) < 2e~ c * E(M) 

(See [1J, Cor A. 14). 
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Since (p n (D n ) converges to D with respect to the Hausdorff distance, 



equation (4.2) guaranties E \T^ nD ) is at least a power of n, for all z E 



Z such that [ [n d + x zA ,n d+1 Zd) E D t n if. Hence, there is a 

\\ / ie{l -d-l} 

constant C2 uniform in z such that: 



E(.M n (*)) > E (£ n (z)j > c 2 nw. 



Hence, both E yC n (z)j and E (M n (z)) are at least powers of n. 

As a consequence, with a probability that is exponentially close to 1, 
uniformly in z such that (n~ d ~+ i Z\...d-i- l ri~ d ~+ i Zdj E D t (~]K, the point z lies 
inside A(n). 

Using a simple union bound on the polynomial number of points z such 

that: ( (rT d ~+ l Zi ) , rC^Zd ) E D e H K, the Borel-Cantelli lemma 

\\ Jie{i-d-i} ) 

states that almost surely, for all but finitely many n, 

D e c K (A(n)) . 

Remark: The odometer function for the iDLA model does not appear 
directly in this proof. However, the function f n ,D n is closely linked to 
it. Hence, there is reason to believe that the normalised odometer for 
the drifted iDLA model also converges towards the difference between the 
obstacle function 7 and its least supercaloric majorant, as suggested by the 
simulation presented in Figure [3] 

Theorem 4.2 (External convergence towards D) Let A{n) be the drifted 
iDLA cluster. 

Then, for every e > ; 

(j) n {A{n)) C D\ 
for all but finitely many n, almost surely. 

First, we remark that the regularity of the solution to the obstacle 
problem is such that the Lebesgue measure of the boundary of D is 0, as 



stated in proposition 3.13 Together with Theorem 4.1, it means that for 
every e > 0, we can find rj > and K such that the number of points z 
such that: 

n^^zA ,n~d^Zd) E D e H K 

J ie{i-d-i} J 

is more than (1 — e)n. 

Using the classical approach of [21j, we argue that we only need control 
over the remaining en points. We will need to show that these particles 
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Figure 3: On the top, the obstacle function j(x, t), and on the bottom, the 
sum of the obstacle function and normalised odometer. The drifted iDLA 
cluster can be glimpsed as the non-coincidence set. 
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do not make the cluster grow beyond distance Ce of its internal limiting 
shape D, with C a constant depending only on the dimension. We will do 
this using a layers argument. 

First, we need to argue that the particles spread tout considerably when 
they travel a macroscopic distance in either time (the direction of the drift) 
or space (other dimensions). In the direction of the drift, this is ensured 
by the local central limit theorem. The case of other dimensions is studied 
in the following lemma. 

For simplicity reasons, we begin by stating the desired result using a 
simple random walk in the (d — 1) space dimensions. 



Lemma 4.3 (First contact point with a height set) Let X n be a (d— 

1)- dimensional simple random walk starting at the origin. Let be the 
hitting time of level k, namely Tk = infjt e N, H-XtH^ > k}. 

Then, given x > 0, there exists a constant C, such that for all x, t G 



(R* + ) , with x > x Q , 



\ [xnd+l] 1 







< Cn d + 1 - 



This result can be explained as follows : the probability P n that tn d + 1 

i 

is the first hitting time of level xn d + l is the product of the probability that 

1 2 

the value xn d + l is taken at time tn^ 1 by the probability, conditioned on 
this fact, that this level is never reached before: 



i = tn d +! 

i+T 



X 



xn d +! 



tnd+i 



When d = 2, a symmetry argument shows that the right-hand term is 
equal to the probability, for a walk starting at the origin and conditioned 
to reach xns at time tns, to stay always strictly positive (except at the 
origin), and never go above 2xn~ d ~+ l . 

Considering only the first condition, our conditional probability is less 
than the probability, for a walk starting at the origin and conditioned to 
reach xna at time tn^ , to stay always strictly positive (except at the origin). 

This is a well-known problem (corresponding to the probability, when 
counting ballots, that a candidate with xnz more votes than the other is 
always first on the tally when the total number of votes is tns), and this 
probability is exactly the following : 



(V < tra§,t' > 0,X V > 



X 2 

tn3 



1 

xn 3 



xn-i 

2~ 

tn 3 



x 



-71 
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When d > 2, we can apply the same argument. Indeed, suppose that 

i 

the random walk reaches level m^ 1 through coordinate i. Then consider 
the projection X 1 of X to this coordinate. It is a lazy random walk, that 
remains in place with probability at each step and does a simple random 
walk otherwise. 

For any k, the probability that ||^||oo stays less that k can be bounded 
above by the probability that \X l \ stays less than k. 

Since X 1 is a lazy random walk, the number of non-zero steps it takes 

2 2 

in time tn^ 1 is asymptotically greater than Atn^ 1 , where A is a suitable 
constant. 

Hence, applying the previous result to X 1 yields : 



i 

X j_ = < 5- 



Using the local central limit theorem, we get: 

P it i = tn^i | < P ( X 2 = xn^Fi ) — n~dii 



Ccc ( 1 1 x 



^ M exp {-^r ]n d+1 



where the constant C is uniform in x > Xq. This concludes the proof of 
lemma 14.31 

Let us now consider our drifted random walk. We consider its projection 
on the space dimensions, with a change of the time parameter such that is 



a simple random walk. Lemma 4.3 then states that the distribution of its 



crossing time for level [xn^+ij is spread out: 

P (t i = \tn^ \) < Cn~^. 

\ [xnd+i] L J J — 

We now argue our drifted random walk shares the same property, in the 
sense that the point at which its trajectory crosses the level xn^ 1 is spread 
out. Let us define r k = inf{t <E N, 3z 6 {1 • • • d - 1}| (S(t)) i | > k}. Then 
there is a constant C uniform in x > Xq > 0, t, such that: 

p (( s ( r 1 - l ))r |, " l,| )- ftr4 



Indeed, the position in time of the crossing point is the number of steps 

i 

taken by S during the time it takes its projection to reach the level xn^ 1 , 
which spreads out its distribution even more. 

We will now bound the maximum distance that one of the en particles 
can reach outside the cluster in the space dimensions. 
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We build the cluster and stop the particles when they exit the discrete 
version of D e C\K, so that we have a number a(n) < en of particles waiting 
to be released on the boundary, and we define A(i) as the cluster when i of 
these particles have been released. It is a well-known property of the iDLA 
models that the law of A(a(n)) does not depend on the order in which they 
are released. We will re-index these walks as (S' l )i e fi... a f n )}- 

Choose t]q > 0, and consider the layers of points 

H k ={ze Z d , [n^ V o\ < d h (z,D) < [n^ VQ \ + l} , 

where dh(z,D) denotes the distance, in Z d_1 , between z and set of points 
that normalise into D in the d — 1 non-drift directions (We ignore the 
distance in the direction of the drift). We will study the expected value of 
the number of particles in each layer T-i k . Le ut us define : 

v H>k (i) =E(\n k nA(i)\ 

When the Z + 1-th random walk S l+1 is launched from the point at which 
it was stuck, consider the event that it adds to the cluster at a point of 
% k . Now if S l+l adds to the cluster at a point of T-L k , it means that S l+1 
crosses %k-i for the first time at a point of Hk-i V\A{f). The distance that 
the random walk has to cross is bigger than r]Qn^+ l , so that using equation 



(4.2) gives a uniform bound the hitting probability of any given point. 1 
simple union bound yields the following recursive relation, for a suitable 
constant C\ : 



m,k(l + !) _ MO < Cin n , h -i{l)n w 

Summing over values /, and remarking that fik-i{0) = 0, yields the 
following equation: 

o(n) — 1 

1=1 

By a simple induction argument, using the inequality Ym=i 1 l k 1 — 
"y , we have the following inequality: 

_d_\ k a(n) k 



yi,H,k{oi{n)) < (Cxn d+ij 



i a{n) d_ x " 

Choosing k = Ken^ 1 , with K suitably big, results in the fact that 
A t K,fc( ct ( n )) is exponentially small. Hence, by the Borel-Cantelli lemma, the 
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normalised cluster contains no points further away from D, in the direction 
of the drift, than r/o + Ke. 

We proceed in the same way to bound the cluster in the direction of 
the drift. 

Choose t] > 0, and consider the layers of points 

V k = {z G Z d , \_rj n^\ < d v (z,D) < [ Vo n^ \ + l} , 

where d v (z,D) denotes the distance, in Z, between z and points that nor- 
malise to D, in the direction of the drift. 

We define //y^(7) as the expected value of the number of particles settled 
in Vk after launching / particles. 

Let us consider a walk S l . Since all points of Vk are at least at distance 

2 

r)on d + 1 from the starting point of S l , the local central limit theorem gives a 
uniform bound for the probability P ra of hitting any point of Vk before all 
the others: for a suitable constant we have: 

P„ < C 2 n^. 

Once more, since a particle needs to cross Vk-i on a point of A(l) in 
order to contribute to A(l + 1) at a point of Vk, we have a recursive relation 
that is similar to that of the previous proof, and for a suited constant C3, 



-1 



»v,k{l + 1) - vy,k( l ) ^ C 3 u Vt k-i(l)n <*+i 

a(n) — 1 



1=1 



By induction, this yields: 



d-l \ Tl^ 

vv,k{Gi{n)) < (a 3 n - 3+T) — 
, , oc\n) _d = i ■ ' 



Choosing k = Ken^, with K suitably big, yields that vy^{a(n)) is 
exponentially small. Hence, by the Borel-Cantelli lemma, the normalised 
cluster contains no points at distance greater than r] + Ke from D. This 



concludes the proof of Theorem 4.2 
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4.3 A more natural class of drifted random walks 



Throughout our proofs so far, we have considered the drifted random walks 
S 3 with the following law: 

P (Si (t + 1) - S> (t) = ±eA = }~ P _,, fori = l-..d-l, and 

2 (a — 1) 



F(S j (t + l)-S j (t) = e d ) = p. 



(4.5) 



A more natural drifted random walk can be defined as a walk that, 
at each time step, performs a simple random walk on Z d with probability 
1 — p, and a step in the direction of the drift with probability p. Such a 
random walk S' has the following law: 

F(S'(t + l)-S'(t) = ±e l ) = X -^- fori = l---d-l, 

¥(S'(t + l)-S'(t) = e d ) = P+^p and 

F(S'(t + l)-S'(t) = -e d ) = 

In order to extend our limiting shape result to the cluster built using 
this random walk, we will see step by step how our proofs need to be 
modified. 

First, we need a different version of the unfair divisible sandpile, in 
which the mass emitted from a point during its toppling is distributed in 
accordance with the new law of our random walk: a fraction ^ is sent to 
each neighbor, and an additional fraction p is sent towards the neighbor in 
the direction of the drift. 

Once again, we normalise our cluster by in non-drift directions, 

2 

and n d + x in the direction of the drift. 

Remark that the new version of the discrete operator we define still 
verifies the minimum principle, so that the final configuration of mass is 
once again defined as the lowest supercaloric majorant of a suitable obstacle 
function, where the term supercaloric is defined with respect to the new 
discrete operator. 

One suitable obstacle function is the following: 

d d ~ l 

= ~J^1 5^ ^ + Zd ~ ng ^ 0, ^' 

i=l 

where g is the discrete Green function for the random walk S'. 

Coupling our new drifted random walk with a random walk on Z 6 ^ 1 with 
a suitable law, we get an convergence of the normalised obstacle function 
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towards the function: 



1 f (3 \ — f lb 2 



1 '(x,t) = t- — \x\ 2 --{-) exp(-/3 

where f3 = 2 (i- P ) > ^ na ^ * s once niore uniform on compacts of ]R d_1 x M* + . 
We also have a suitable bound on the error. 

The next step of the proof is to adapt the continuous operator .ft, which 
becomes: 



Then, following the proof of lemma [379| , we can prove that the normalised 
version of the least supercaloric majorant of 7 n converges to the least su- 
percaloric majorant of 7', where the supercaloric functions are defined as 



in section 3.3.2 with respect to the new heat operator: 



* /_ 2d A/ P df 

This convergence is once again a little technical to prove, but it relies 
only on the precise estimates of the convergence of the Green function, 
and the fact that our discrete operator is a good approximation of the 
continuous one for functions with sufficient regularity. Remark that the 
terms corresponding to a simple random walk in the direction of the drift 
become negligible because of our normalisation. 



Using the same arguments as in the proof of Theorem 3.1 we can thus 
prove that the normalised cluster for the new unfair divisible sandpile 
model, when intersected with a compact, converges towards the limiting 
shape D D K, where D is defined as the non-coincidence set of 7' and its 
least supercaloric majorant. 

The proof of the extension of the convergence result to the iDLA cluster 



follows that of Theorem 4.1 almost verbatim, and the majoration of the 



exterior error can be derived from the proof of Theorem |4.2[ using similar 
estimates on the probability for a random walk of hitting any given point 
in a set, after covering a macroscopic distance in any direction. 



5 Properties of the cluster 

We know that our cluster converges to a limiting shape D defined in terms 
of an obstacle function, however we are interested in the properties of this 
limiting shape. To natural questions arise, the first being its boundedness. 
Indeed, our normalisation only truly captures the behaviour of the cluster if 
we can prove that the limiting shape under this normalisation is bounded. 
Moreover, we know that the limiting shape is a true heat ball, and the 
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existence of a bounded true heat ball is a question of interest in PDE 
theory. 

The second natural question concerning the shape D is that of its uni- 
versality: when we change the parameters of the model, we will show that 
the limiting shape is only modified through a simple variable change. 



5.1 Bounds on the cluster 

In this section, we work once again with our drifted random walk S, that 
performs at each step a simple random in Z d with probability 1 — p, and 
steps in the direction of the drift with probability p. Using probabilistic 
arguments yields bounds on the iDLA cluster, which can in turn be trans- 
ferred to the continuous shape D. Hence, our limiting shape D will be a 
bounded true heat ball. 



5.1.1 Horizontal bound 

Lemma 5.1 (Non-drift direction bound on the cluster) The iDLA clus- 
ter A(n), normalised by n 3 ? 1 , is bounded in all "non-drift" directions. 



Proof: We proceed as in the proof of Theorem |4.2 and we consider the 



intersection of the entire cluster with strips in successively non-drift and 
drift directions. 

Choose xq > 0, and define ko = XqU 1 ^ , and Tk = {(x, t) G IR d_1 x 
R, | \x\\oo = ko + k} In this proof, we study the expected value of the number 
of particles in each cylindrical domain I\.. Define /ifc(Z) = E (A(l) R Tfc). 

When the / + 1-th random walk X l+1 is launched from the origin, con- 
sider the event that it adds to the cluster at a point of IV Now if X l+1 
adds to the cluster at a point of T^, it means that X l+1 crosses T^-i for 



the first time at a point of Tk PI ^4(0- Lemma 4.2 gives a uniform bound 
on this hitting probability, so that we have the following recursive relation, 
for a suitable constant C\ : 

fi k (l + 1) - MO < C 1 /X fc _i(Z)n"3¥i 

Summing over values I, and remarking that /i J t_ 1 (0) = 0, yields the 
following equation: 

n-l 

d 



l=i 
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By a simple induction argument, using the inequality Ym=i ^ 1 — IT' 
we have the following inequality: 

< (=<*»-*)' 

Choosing = xin^+i, with xi suitably big, results in the fact that /J, k (n) 
is exponentially small. Hence, by the Borel-Cantelli lemma, the normalised 
cluster is almost surely asymptotically bounded by x + x±. 



5.1.2 Vertical bound 

Lemma 5.2 (Drift-direction bound on the cluster) The iDLA clus- 

2 

ter normalised by n d + x in the direction of the drift is bounded in the direc- 
tion of the drift. 

Proof: The Local Central Limit theorem gives a uniform bound for the 
hitting probability of a vertical strip at a particular point: for a suitable 
constant C 2 , we have: 



/ i \ d-i 

[X 2 = in^ 1 < C 2 n d + 1 



We will follow a similar strategy as in the previous proof. 

2 

Choose t > 0, and define k = t n d +! and the vertical strips V k = 
{(x, t){(x, t) G IR^" 1 x R, t = k + k}. 

We will study the number of particles in each vertical strip V k . Let 
z/fc(7) = E (A(Z) n Ck +k) denote its expected value. 

Since a particle needs to cross V k _i on a point of A(l) in order to con- 
tribute to A(l + 1) at a point of V fc , we have a recursive relation that is 
similar to that of the previous proof, and for a suited constant C3, 

v k (l + 1) - v k (l) < C 3 i/ fc -i(Z)n"a^ 

n-l 



By induction, this yields: 

/ d-i \ 1 
v k {n) < (C 3 n d+ij 



z=i 



< 



k\ 

k 



Choosing k = tin^ 1 , with t\ suitably big, yields the summability of 
z/ fc (n). Hence, by the Borel-Cantelli lemma, the normalised cluster is almost 
surely asymptotically bounded by to + t±. 
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5.2 Rescaling 



In this section, we will show that the various limiting shapes can be written 
in terms of one another by rescaling differently in the direction of the drift 
and in the other directions. 

Lemma 5.3 Let pi and p 2 be two drift parameters in (0, oo), and let D\ 
(respectively D 2 ) be the limiting shape of the unfair divisible sandpile model 
run with parameter pi (respectively p 2 ) ■ 

Then D 2 is the image of Di by a change of variables x — > jix,t — > At. 

First, we introduce the additional parameter k G which measures 
the quantity of mass sent from the origin. We will call unfair divisible model 
with mass k the unfair divisible sandpile model run with initial mass kn 
at the origin. Remark that the limiting shape of this model is of course 
obtained from that of the original model by a rescaling (with different 
coefficients in the drift and non-drift directions). Moreover, the limiting 
shape of the model can still be obtained as the solution to our parabolic 
obstacle problem, the only difference being an adjusted coefficient k in the 
last term of the obstacle function: 



/ x , ,i2 k f (3\ 2 * ( \\x\ 

7 *0M) = t-\\x\\ exp^-/3 — 

We are now ready to compare D\ and D 2 . Since the shape D stems 
from the obstacle 7, we only need to show that 72, the obstacle function 
for parameter p 2 can be related to 7^, the obstacle function for drift p\ 
and mass k. Consider the two following functions: 

1 / x x A w 2 ,, M2 k ( ' fixX^ ( n X \\x\ ' 

72OM) = *-|MI 2 -^(|) 2 ex P (-&M 

with C a suitable constant such that 

AI^-g||WN=-l. 

We see that the constant f3i can be changed to (5 2 = provided that 

k is such that the last term of ^71,* is equal to that of 72, that is to say, 

fci 1 



Cp^V p 2 
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Since ^jx^fxx , At) and 72(3?, t) now only differ by a caloric function, they 
give rise to the same limiting odometer, hence to the same limiting shape. 
We conclude using the first part of the proof to recover that D\ and D2 
are indeed images of one another by a transformation of the required form 
x —7- fix, t — > \t. 



6 Conclusion 

At this point, we have proved the convergence of our two models towards 
a limiting shape S that solves the following PDE problem: given <p a C°° 
function of time and space such that 

Acp + p— = 0. 



2(d- 1) r ' dt 
Then the following mean value property holds: 

0(MMM) = 1^10(0). 



s 



We have also proved, using probabilistic estimates on random walks, 
that the iDLA cluster is bounded. Moreover, the regularity of the problem 
enables us to transpose these bounds to S, since S is sufficiently defined in 



terms of a Hausdorff limit. This concludes the proof of Theorem 1.1 
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